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Abstract 

We study the full contributions at the leading order(LO) and QCD next-to-leading or- 
der(NLO) to the cross section of the e + e~ — > W + W~~bb process in the standard model(SM) at 
the ILC. In dealing the resonance problem we adopted the complex mass scheme in both tree- 
level and one-loop level perturbative calculations. Our numerical results show that the K-factor 
varies from 1.501 to 0.847 when </s goes up from 360 GeV to 1.5 TeV. We investigate the 
dependence of the LO and QCD NLO corrected cross sections of process e + e~ — > W + W~bb on 
colliding energy y/s and Higgs-boson mass. We also present the results of the LO and QCD NLO 
corrected distributions of the transverse momenta of final particles, and the invariant masses of 
Wb-, bb- and WW-pait. 



PACS: 13.66.Jn, 14.65.Ha, 14.80.Bn, 12.38.Bx 
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I. Introduction 



The Higgs boson, which gives masses to the weak vector bosons and fermions, plays an important 
role in the standard model(SM). Unfortunately, it has not been directly detected yet in experiments. 
Searching for Higgs boson within the standard model(SM) and study the phenomenology concerning 
Higgs properties are the important tasks at the present and upcoming high energy colliders. LEP II 
experiments have provided the lower limit on the SM Higgs mass as 114.4 GeV at the 95% confidence 
level, which is extracted from the results of searches for e + e~ — > Z°H° production [H [2J. While the 
indirect evidences of the SM Higgs mass through electroweak precision measurements indicate the 
95% C.L. upper bound as mu < 182 GeV , when the lower limit on mu is used in determination of 
this upper limit [2J. On the other hand, the heavy top-quark practically plays a central and crucial 
role in probing the electroweak symmetry breaking as well as the flavor problem in all the extended 
models beyond the SM which address the hierarchy problem. Recently, a new datum of top-quark 
mass has been already presented by the CDF and DO experiments at Fermilab, and the preliminary 
world average mass of the top-quark is known as mt = 172.5 ± 1.3(stat) ± 1.9(syst) GeV, which 
corresponds to a 20% precision improvement relative to the previous combination [3J. 

The future International Linear Collider (ILC) is proposed by the particle physics community 
with the entire colliding energy in the range of 200 GeV < y/~s < 500 GeV and an integrated 
luminosity of around 500 in four years. The machine should be upgradeable to y/s ~ 1 TeV 

with an integrated luminosity of 1 (ab)^ 1 in three years [1]. Most of the main physics topics within 
the SM or its beyond at TeV energy scale can be explored at such a machine. Emphasis is given 
to the study of top-quark physics, electroweak physics in the SM, and the measurements in the 
extended SM, such as supersymmetry. 

Compared with the hadron colliders, such as the Tevatron and the CERN Large Hadron Col- 
lider (LHC), the ILC can produce top and Higgs boson signal events more easily resolved from 
backgrounds. Therefore, the ILC is an ideal facility to study top and Higgs physics with much 
more precise measurement for their parameters. At the ILC we can also carry out the study of 
gauge boson interactions, and the delicate cancellations which are related to the gauge structure of 
the theory and essential to preserve unitarity. Furthermore, the ILC experiment might be able to 
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explore the signature of the new physics, if the SM is really only an effective theory at low energy. 

At the ILC, detecting the top-quark pair production process e + e~ — > ti is a good way to study 
the top-quark properties, and the associated Higgs production with Z boson e + e~ is 
one of the cleanest signature in discovering Higgs boson if the the b-quark trigger system has high 
performances except vertex detectors [5]. The former process will be followed by the subsequential 
decay through ti — ► W + W _ 66[B], while the later process goes via H°Z° — ► W + W~bb through 
decays H° — ► W + W~ and Z° — ► bb if the Higgs boson mass is larger than 2mw[7J. Therefore, 
the signature of e + e~ — ► W + W~bb at the ILC serves as non-resonant background to both top- 
quark pair production and associated production of Higgs boson with Z° boson. We can see 
that it is crucial to separate the top and Higgs signatures from the other W + W~bb production 
backgrounds in ILC experimental data analyzing. In the precise measurements of the signals of 
both the ti pair and H°Z° associated production processes, the relevant irreducible background 
from e + e~ — > W + W~bb should be carefully investigated. 

In Refs. [51 151 H?l [TU 1 \TT \ fT2"] the NLO electroweak and QCD corrections to the process e + e~ — ► ti 
and decay t — * W + b have been already extensively studied. And the non-relativistic effect near the 
threshold of tt production is also studied carefully in Ref.[13j, which can not be reliably described 
with fixed QCD orders in perturbative theory. The Higgs-strahlung Bjorken process e + e" 
was investigated in Ref. [14] . and the process e + e~ — > tt — > W + W~bb — > 6f with six fermion final 
states after W pair decays has been also calculated at the lowest order in Ref. [15] . The evaluation 
of the e + e~ — ► W + W~bb process with finite width method at the tree-level is also presented in Ref. 
[71 [in] . All those studies indicate that the precise investigations of the characteristics of top-quark 
and the Higgs-boson are significant for the future e + e~ ILC experiments. 

In this paper we present the calculations of the cross section of the process e + e~ — > W + W~bb 
at the leading order(LO) and its QCD next-to-leading order(NLO) (0(a s )) corrections. The paper 
is organized as follows: In the following section we present the analytical calculations for process 
e + e~ — > W + W~bb at the LO and QCD NLO. The verifications of the correctness of our calculations 
are declared in section III. The numerical results and discussions are given in section IV. In the 
last section we give a short summary. 
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II. Calculations 

The calculations for the process e + e~ — > W + W~bb are carried out in 't Hooft-Feynman gauge. 
In the QCD NLO calculations, we use the dimensional regularization(DR) method to isolate the 
ultraviolet (UV) and infrared(IR) singularities. In order to preserve gauge invariance, we adopt the 
approach of the complex mass scheme to deal with the unstable particles in the calculations for the 
tree-level cross section and QCD NLO radiative correction [171 \TE\ . The on-mass-shell(OS) scheme 
is used to renormalize the masses and fields of related bosons and fermions. The FeynArts3.2 
package [TU] is adopted to generate Feynman diagrams and convert them to corresponding ampli- 
tudes. The amplitude calculations are mainly implemented by applying FormCalc4.1 programs[20j. 
The formula for calculating the IR divergent integrals with complex internal masses in DR scheme 
are obtained by analytically extending the expressions in Ref.[2T] to the complex plane. The numer- 
ical evaluations of IR safe one-point, two-point, three-point and four-point integrals with internal 
complex masses, are implemented by using the expressions analytically continued to complex plane 
from those presented in Refs.[22J [23]. And the 5-point scalar integral can be expressed in terms 
of multiple scalar four-point integrals |24j. The subroutines for one-loop integrals with complex 
masses are coded based on the LoopTools2.1[20| package which comes from FF library [2"5]. The 
2 — ► 4 phase-space integration routine[27] is created based on the 2to3.F program in FormCalc4.1 
package. The five-body phase-space integration for hard gluon radiation process e + e~ — ► W + W~bb 
g is accomplished by using CompHEP-4.4p3 program |26j. 

Now we present the analytically calculations of the tree-level cross section for e + e~ — ► W + W~bb 
and its QCD NLO radiative corrections. The notations for the process are defined as 

e + (pi) + e-(p 2 ) - W + ( P3 ) + W-(p A ) + 6(p 5 ) + b(p 6 ), (2.1) 

where Pi (i = 1 — 6) label the four- momenta of incoming e + , e~ and outgoing final particles, 
respectively. There are 64 generic tree- level diagrams for the process e + e~ — ► W + W~bb presented 
in Fig(TJ where internal wavy- line represents 7, Z°, or and internal dash- line represents a 
Higgs-boson H° or a Goldstone G (G =t ). We can easily find that in FigfT] there includes the tree- 
level diagrams for the processes e + e~ — ► t*t* — > W + W~bb and e + e~ — ► H°*Z°* — ► W + W~bb. 
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Figure 1: The generic tree-level Feynman diagrams for the e + e — ► W + W bb process. Internal 
wavy-line represents 7-, Z -, or VF^-propagator. Internal dash-line represents a Higgs boson H° 
or a Goldstone G°(G ± ). 



The differential cross section for the process e + e — ► W + W bb at the tree-level is obtained by the 
tree-level is obtained by 

(2vr) 4 



d<7 



tree 



^2\M t ree\ 2 d<f> A , 



(2.2) 



where d^n is the four-body phase space element given by 

/ 6 \ 6 



d$ 4 = 5^ [ P i+P2- Z zpi n 



d 3 pi 



(2tt) 3 2^ ' 



(2.3) 



\ i=3 / i=3 

The summation in Eq. (|2,2p is taken over the spins and colors of final states, and the bar over 
the summation recalls averaging over initial spin states. In the calculation, the internal Z° and 
Higgs boson can be real, and the top-quark propagator can also be resonance when y/s > 2mt- 
To deal with these resonant singularities, we use the so-called complex mass scheme(CMS) in our 
perturbative calculations |1 7\ [T8], The complex masses of W-, Z-, H-boson and top-quark are defined 
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as 



(X = W,Z, H,t). 



(2.4) 



In the CMS approach the complex masses for all related unstable particles should be taken ev- 
erywhere in both tree-level and one-loop level calculations. Then the gauge invariance can be 
conserved and singularity poles of propagators are avoided. 

In calculating the complete QCD NLO corrections, we should consider the contributions of 
30 self-energy diagrams, 94 triangle diagrams, 17 box diagrams and 6 pentagon diagrams. As a 
representative selection, we present the pentagon Feynman diagrams of the e + e~ — > W + W~bb 
process in FigJ2J We adopt the Eqs.(4.26) and (4.27) in Ref.[TB] for the renormalized QCD self- 
energy and counter-terms of top-quark with complex mass neglecting terms of 0{a 2 s ) by using 
OS-scheme. 

II'; 
W ' 




Figure 2: The pentagon Feynman diagrams for the e + e — * W + W bb process. 



There exist both ultraviolet (UV) divergency and infrared(IR) soft singularity in the contribu- 
tions of the QCD one-loop diagrams for e + e~ — ► W + W~bb process, but no collinear IR singularity 
due to the massive top- and bottom-quark. After doing the renormalization procedure, the UV 
singularity is vanished. 

To cancel the IR soft divergency appeared in the virtual correction, we should consider the 
contribution of the real gluon emission process e + e~ — > W + W~bbg. We denote the real gluon 
emission process as 

e+(pi) + e"(p2) -» W + ( P3 ) + W-(pi) + 6(pb) + Hpg) + g(pr)- (2-5) 
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To calculate the contribution of this process, we introduce an arbitrary small soft cutoff S s to sepa- 
rate its 5-body phase-space into two regions[2B3, i.e., soft (£7 < 5 s ^fs/2) and hard(i?7 > 5 s ^fs/2) re- 
gions. After adopting the soft gluon approximation, the expression of a so ft for e + e~ — ► W + W~bbg 
process with soft gluon has the form as 

dOsoft = Cf^556 da t ree, (2.6) 

where Cp = 4/3 and g§Q are defined as: 

1 2(p 5 • p e ) 



4(m) ln(a) + 2 



e A 1 /2( S56 , m 2 )m 2^ 



A 1 /2( S56 , m 2 ;m 2) 

x [ln 2 (a)+4Li 2 (l-^)] - - lncr + 0(e). (2.7) 



P 

In above equation, \{s^,m 2 ,m 2 ) is the kinematical function defined by: 



X(x, y, z) = x 2 + y 2 + z 2 - 2xy - 2yz - 2zx, (2.8) 



AE = E 7 = 5 s ^/s/2, s 56 = (p 5 + j> 6 ) 2 and 



(2.c 



«56 1 + P 

Our created 2^4 phase space integration routine[27j, is adopted in the tree-level and one-loop 
level calculations for e + e _ — > W + W~bb process. The IR singularity part of the soft gluon emission 
process e + e~ — * W + W~bb (g) can be exactly cancelled by the IR singularity induced by the one- 
loop virtual gluon correction. We apply CompHEP-4.4p3 program [26J to implement the phase space 
integration of the hard gluon emission process e + e~ — * W + W~bb +g. Finally, we get the finite 
total cross section including complete NLO QCD corrections for the process e + e~ — ► W + W~bb by 
summing up all the contribution parts, 

&NLO = &tree + ^virtual + & soft + Chard- (2-10) 

III. Checks 

We have performed the following checks to prove the reliability of our calculation: 
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m t (GeV) 


a LO (fb)(Rel[l6\) 


olo (/b)(Comphep) 


<?Lo{fb) (ours) 


150 
180 
200 


663.11 
576.26 
497.63 


663.03(1±0.05%) 
576.19(1±0.04%) 
497.58(1±0.04%) 


663.19(1±0.05%) 
576.52(1±0.04%) 
497.68(1±0.04%) 



Tabic 1: The comparison of the numerical results of LO cross section neglecting the 
diagrams with Higgs-boson interchanging by using CompHEP-4.4p3 system, our in-house 
2 — ► 4 phase-space integration routine with the corresponding selected results presented 
in Ref.pl] when ^ = 500 GeV. 

• The LO cross section for the process e + e~ — > W + W~bb is calculated in the conditions of 
taking y/s = 500 GeV and neglecting the contribution of the diagrams with internal Higgs- 
boson exchange which are taken in Ref.|16j. The numerical results of the LO cross section 
for the process e + e~ — > W + W~bb are listed in Table [TJ There our results are obtained by 
using both CompHEP-4.4p3 program and our created 2^4 phase-space integration routine, 
and compared with the corresponding ones presented in Ref.|16j. We can see there is a good 
agreement between ours and those presented in Ref.[16]. The in-house 2 — ► 4 phase-space 
integration routine was also once verified in our previous work[27j. 

• We use our created codes for numerical evaluation of the one-loop integrals with complex 
internal masses. The comparisons are made between the results and those obtained by doing 
directly the integration of Feynman-parameter. There exists a good agreement. The results 
from both calculations for scalar two-, three-, four-point integrals are coincident with each 
other at least up to six digits, respectively. 

• The exact cancellations of UV- and IR-divergencies are verified both analytically and numer- 
ically in our calculation. 

• The independence of the total cross section including the NLO QCD corrections on the soft 
cutoff S s (= 2 AE/yfs) is confirmed numerically. Our calculation shows the errors of the 
independence are less than 0.6% in the 5 S region of [10 -4 , 5 x 10~ 2 ]. In further numerical 
calculation we fix 5 S = 10~ 3 . 

• In the following section, we shall clarify other verifications. 
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IV. Numerical results and discussion 



In our numerical calculation we take the following input parameters 

a{m z )- 1 = 127.918, a s (m|) = 0.1176, m z = 91.1876 GeV, 

m w = 80.403 GeV, T z = 2.495 GeV, T w = 2.141 GeV, (4.1) 

m t = 172.5 GeV, m b = 4.7 GeV, m e = 0.5109991 MeV 

Due to the application of the CMS approach, we use the complex weak mixing angle defined as 

£ = 1-4 = 4. (4.2) 

In our LO and NLO numerical calculations we set the QCD renormalization scale \x as [i = mw+mb, 
and take the strong coupling a s (/i 2 ) = 0.11885, which is obtained by using the formula at three-loop 



level {MS scheme) with the five active flavors [29 . 

Since the widths of top-quark and Higgs boson haven't been well provided or measured exper- 
imentally by now, we use their theoretical results from perturbative calculations. Considering the 
fact that top-quark mass is above m\v + rrif,, and V#, ~ 1, the decay of top-quark is dominated by 
undergoing two-body decay t — ► W + b, and the total decay width of top-quark is approximately 
equal to the decay width of t — > W + b. Neglecting terms of order m^/mf, a 2 s and (a s /7r)M^y/m 2 , 
the width predicted in the SM is [31] : 



2 / I ' 2 

171$ J \ m,t 



(4.3) 



16m^(l - mfr/mf.) 

The reasonable physical decay width of Higgs boson is obtained by employing the program Hdecav[32]. 
where the partial decay width F(H° — > qq) is calculated including O(af) QCD radiative cor- 
rections. Then we obtain r t = 1.3745 GeV, T H (m H = 120 GeV) = 0.3692 x 10~ 2 GeV and 
T H {m H = 180 GeV) = 0.6286 GeV. 

The numerical results of the LO, QCD NLO corrected cross sections and the corresponding K- 
factor(if = a ^° ) for the process e + e~ — * W + W~bb are plotted in Figsl^a) and (b) respectively, 
when ran = 120 GeV. As indicated in Figj3|a), both curves for the cross sections at the LO and 
NLO increase quickly in the y/s region of [350 GeV, 400 GeV] and decrease when ^fs > 430 GeV. 
Figj3jb) shows that the corresponding K- factor decreases slowly from 1.501 to 0.847 as y/s running 
from 360 GeV to 1.5 TeV. The large positive peak near the ti threshold in Figj3]^b) is due to a 
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Figure 3: (a) The LO and NLO QCD corrected cross sections for the process e + e~ — ► W + W~bb as 
the functions of c.m.s. colliding energy(- v /s) with m# = 120 GeV, (b) the corresponding K-factor 
versus y/s. 



y/s{GeV) 


&tree(fb) 


VNLo(fb) 


K-factor 


500 

1000 

1500 


602.57(1±0.05%) 
182.24(1±0.04%) 
82.73(1±0.04%) 


575.5(1± 0.38 %) 
156.7(1± 0.38 %) 
70.1(1± 0.37 %) 


0.955(4) 
0.860(4) 
0.847(4) 



Table 2: The LO and NLO QCD corrected cross sections, K-factors for e + e — > W + W bb 
process with m H = 120 GeV and y/s = 500 GeV , 1000 GeF, 1500 GeV , respectively. 

Coulomb singularity effect coming from the instantaneous gluon exchange between heavy quarks 
which has a small spatial momentum. In Table [2] we list the values of <Tt ree , &nlo an d K-factor 
at some typical y/s points, which are read out from Figsj3ja-b). Since the QCD correction to the 
e + e~ — * W + W~bb process with high colliding energy can be approximately decomposed into the 
QCD correction to the tt production plus the corrections to the t(t) — ► W + b(W~b) decays when 
mjj < 2mwi we make following verification to check our results. We evaluate the QCD correction 
to e + e~ — > W + W~bb process by combining the QCD corrections to e + e~ — > ti production and 
t(t) — ► W + b(W~b) decays together, and get the K-factors to process e + e~ — > W + W~bb as 0.8562(1) 
for y/s = 1 TeV and 0.8433(1) for y/s = 1.5 TeV, which are coincident with the corresponding 
ones in Table [2] in error ranges. 

In Fig|U(a) we present the plot of the LO and QCD NLO corrected cross sections as the functions 
of Higgs-boson mass, with y/s = 500 GeV and run running form 60 GeV to 200 GeV . We find from 
FiglU(a) that the LO and QCD NLO corrected cross sections are non-sensitive to the Higgs-boson 
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(b) 



100 120 140 160 180 200 



Figure 4: (a) The LO and QCD NLO corrected cross sections for the process e + e~ — ► W + W~bb as 
the functions of Higgs mass(m# ) with = 500 GeV . (b) The corresponding relative QCD NLO 
corrections versus ran- 

mass except in the vicinity of m# ~ 2mw, from there the Higgs mass becomes larger than 2mw, and 
H -, Z°-boson are simultaneously resonances. We can see also from the figures that the contribution 
via e + e~ — > t*t* — > W + W~bb channel is much larger than that from e + e~ — > H°*z°* — > W + W~bb 
as concluded in Ref.[!l6]. Fig|4jb) shows the corresponding K-factor has the values around 0.956 in 
the range of m H e [60 GeV, 200 GeV]. 

Due to the CP-conservation, the distributions of transverse momenta of VF _ -boson and b quark 
should be the same as those of p>p + and p h T , respectively. We only present the distributions 
of the pj? + and p h T with run = 120 GeV and \fs = 500 GeV in Figsj^a) and (b). In these 
two figures we can see that the QCD NLO corrections suppress the LO differential cross sections 
d&Lo/dpjf + and da lo I dp\ ■ They also show that the differential cross sections of d<JLO,NLO I 'dp^ + 
and d<JLo,NLo/dpj> have their maximal values at about pYp + ~ 70 GeV and p\ ~ 30 GeV respec- 
tively. We see that the line shapes of the differential cross sections in these two figures are mainly 
determined by the contributions of the of e + e~ — > t*t* — » W + W~bb production. 

We plot the invariant mass distributions of (VF + 6)-pair, denoted as M( W + b ^, at the LO and 



QCD NLO in FigE]with m H = 120 GeV and ^/s = 500 GeV. The distribution of M. 



(W~b) 



should 



be the same as that of (W + b)-paii due to the CP-conservation. We can see from the figure that 
most of the events are concentrated around a peak located at the position of M^ w + b ^ ~ nit. That 
demonstrates again the main contribution to the cross section of the process e + e~ — > W + W~bb 



11 




p™(GeV) pJ(GeV) 



Figure 5: The distributions of the transverse momenta of W + and bottom-quark for the e + e — > 
W+W'bb process at the LO and QCD NLO with y/s = 500 GeV and m H = 120 GeV. (a) for 
W + , (b) for bottom-quark. 




140 160 180 200 220 240 



Figure 6: The distributions of the invariant mass of (W + 6)-pair (or (W~b)-p&ii) at the LO and 
QCD NLO with m H = 120 GeV and yfl = 500 GeV. 

with high colliding energy, is from top-pair production channel e + e~ — > tt and followed by the 
decay of t(t) — > Here we can see that the QCD NLO correction slightly suppresses 

the LO differential cross section dcjLo/d^(w+b)- 

As we know, if Higgs boson has a mass larger than 2mw, the e + e~ — ► H°*z°* — ► W + W~bb 
channel will certainly slightly increase both the LO and QCD NLO corrected cross sections for 
W + W~bb production due to the Higgs-boson resonant effect as shown in FigUJa). It will bring 
a spike on the distribution of the (T / l^ + Vl^ _ )-pair invariant mass at the position of M(ww\ = ma- 
Analogously, the associated real Z°-boson produced via e + e~ — > H°Z° will induce a spike on the 
distribution of the invariant mass Mn^\ at the position of Mn^ = mz- In FigsfTJ^a) and (b) we 
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m H = 180 GeV 









Figure 7: The distributions of the invariant masses of (66)-pair and (WW )-pair with mu = 180 GeV 
and y/s = 500 GeV. (a) is for (VFVF)-pair, (b) is for (bb)-pair. 

show the the distributions of the WW- and 66-pair invariant masses with yfs = 500 GeV and 
rriH = 180 GeV , respectively. We can see spikes around the vicinities of M( WW ^ ~ mu ~ 180 GeV 
and Mn^\ ~ mz ~ 90 GeV in FiglT^a) and FigJT^b) respectively, which may be used to distinguish 
the e + e~ — > H°Z° — > W + W~bb signature from the corresponding irreducible background e + e~ — > 
W + W~bb . It shows also the QCD NLO correction obviously modifies the LO differential cross 
sections of daio/dM^ ww ^ and daLo/dM^y 

Theoretically, the NLO QCD correction to the process e + e~ — > Z°H° — > W + W~bb with real Z°- 
and Higgs-boson as intermediate particles, should be determined only by the NLO QCD corrections 
to the Z° — > bb decay. As a check to verify our calculations, we also calculate the correction to 
the decay Z° — > bb with yfs = 500 GeV and mu = 180 GeV > 2mw and get the K-factors 
for e + e~ — > Z°*H°* — > W + W~bb process being 1.0466(1), which is coincident with the result by 
calculating the e + e~ — > Z°*H°* — > W + W~bb process with full QCD NLO diagrams, where the 
K-factors are 1.046(2). 



V. Summary 

In this paper we calculate the complete one-loop QCD corrections in the SM to the process e + e~ — > 
W + W~bb at the ILC. We study the dependence of the LO and QCD NLO corrected cross sections 
of process e + e~ — » W + W~bb on colliding energy y/s and Higgs-boson mass. We investigate the 
LO and QCD NLO corrected distributions of the transverse momenta of final particles and the 
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LO and QCD NLO corrected differential cross sections of invariant masses of Wb-, bb- and WW- 
pair. It shows that NLO QCD correction obviously modifies the LO cross section of the process 
e + e _ — ► W + W~bb , and when the colliding energy ^/s goes up from 360 GeV to 1.5 TeV, the K- 
factor varies from 1.501 to 0.847. The numerical results show that if m# > 2m\y, the resonant effect 
of £f°-boson appearing in the e + e~ — ► H°Z° — ► W + W~bb channel will induce a little enhancement 
to the LO and QCD NLO corrected cross sections for e + e~ — ► W + W~bb process. We find that it 
may be possible to select the e + e~ — > H°Z° — > W + W~bb events from the corresponding irreducible 
background e + e~ — ► W + W~bb which is dominantly produced by the e + e~ — ► tt — ► W + W~bb 
channel by analyzing the invariant masses of final WW- and 66-pair. 
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